General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



X-551-70-313 


NASA TM 


X,fcS32.1 


AN ANALYTICAL DEVELOPMENT OF THE 
RELATIVE MOTION OF TWO CLOSE 
SATELLITES OF AN OBLATE PLANET 


RICHARD W. BARBIERI 


| \r < VI ' 


AUGUST 1970 


0 ° r>l 

•/ 

. r y 

. „ ■} 


GODDARD SPACE FLIGHT CENTER 


N70-39 


GRglgEg, MARYLAND 


(ACCESSION NUMBER) 


YA 


X-55 1-70-313 


A.N ANALYTICAL DEVELOPMENT OF THE RELATIVE' MOTION O^ 
TWO CLOSE SATELLITES OF AN OBLATE PLANET 


Richard W. Barbieri 


AUGUST 1970 


GODDARD SPACE FLIGHT CENTER 
GREENBELT, MARYLAND 


TABLE OF CONTENTS 

Page 

ABSTRACT v 

GLOSSARY OF SYMBOLS vi 

INTRODUCTION 1 

FORMULATION OF THE EQUATIONS OF RELATIVE MOTION 2 

SOLUTION OF EQUATIONS GOVERNING PLANAR MOTION 8 

CONCLUSION AND COMMENTS 14 

REFERENCES. 15 

APPENDIX 1 16 

APPENDIX II 18 

APPENDIX III 20 

APPENDIX IV . 23 




p % c j c 

k£CBD’.HO 


t-AOT 


f.\LKA£^' 


AN ANALYTICAL DEVELOPMENT OF THE RELATIVE MOTION OF 
TWO CLOSE SATELLITES OF AN OBLATE PLANET 


ABSTRACT 

The purpose of this study is to describe the relative planar motion of two 
satellites which are close (■> 5 km) in comparison to their respective distances 
(» 20,000 km) from the center of an oblate central body. This research was under- 
taken to support current studies in orbiting long baseline interferometry systems 
like, for example, the proposed Radio Astronomy Explorer C & D satellite. The 
equations which govern the motion of one of the satellites with respect to the other 
were derived from a Lagrangian formulation and neglect the mutual attraction of 
the two satellites. 

The equations have been expanded to first order in eccentricity of the orbit 
of the reference satellite and also to first order in a small parameter which is a 
function of the oblateness of the central body. The expansion to first order of 
eccentricity is compatible with the intended small eccentricity (< 0.01) of the 
proposed RAE C & D. 

The results to first order eccentricity, neglecting oblateness, show that the 
motion is almost periodic. When the first order oblateness is included a secular 
term appears and the almost periodic behavior becomes superimposed upon a 
secular drift. Certain small divisors appear when the oblateness is included; in 
particular the rate of change of the argument of perigee appears in the form of 
linear combinations with multiples of the mean motion and the modified mean 
motion. 
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Glossary of Symbols 

radius vector of the M satellite with respect to the center of mass of 
the central body 

radius vector of the 0' satellite with respect to the center of mass of 
the central body 

position vector of M with respect to 0' 
the semi-major axis of the orbit of 0' 
the inclination of the orbit of 0' 
the eccentricity of the orbit of 0' 
the longitude of the ascending node of the orbit of O' 
the argument of perigee of 0' 
the true anomaly of 0' 
the sum ( w + v )• 
the angular velocity of 0' 
the kinetic energy of the M satellite 
the potential energy of the M satellite 
the Lagrangian, T-U. 

the velocity components of 0' in the rotating coordinate system. 

the force components, in the rotating coordinate system, acting on 
theM satellite. 

the gravitational constant of the central body 
the equatorial radium of the central body 

the first non zero coefficient in the spherical harmonic expansion of 
the gravity field of the central body 
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n ■ ffj. a -3 
n * n + Qcos t 
M = nt 

w Q the secular rate of change of the argument of perigee 
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AN ANALYTICAL DEVELOPMENT OF THE RELATIVE MOTION OF 
TWO CLOSE SATELLITES OF AN OBLATE PLANET 


INTRODUCTION 

The solution of the problem ol relative motion of two spacecraft has usually 
been constructed in the past with rendezvous and docking applications in mind. 

Since such applications were usually performed within several revolutions about 
the central body it was unnecessary to consider the oblateness of the central body. 
In this report the non-spherical nature of the central body is taken into account 
because problems of a different nature than those of rendezvous and docking have 
received considerable attention over the last few years. In particular, the feasi- 
bility of using orbiting long baseline interferometry systems has been extenshely 
studied. Such a system will be orbiting the central body for a few years at least 
and it would therefore be advantageous to be able to approximate the accumulative 
effects of the oblateness on the motion. Fortunately the intended eccentricity of 
such interferometry systems is small, less than 0.01, and therefore an expansion 
to the first order of eccentricity will be presented. 

The technique used here is rather well known among those working with 
perturbation theory although its application to this particular problem is novel. 

To be explicit, a double expansion has been used to isolate the first order effects 
due to eccentricity from the first order effects due to the presence of an oblate 
central body. Furthermore, even though a device first used by Poincare, that of 
expanding the frequency in terms of a small parameters, hasb» r, n adopted in this 
problem to eliminate mixed secular terms, a small purely secular drift cannot 
be removed from this first order theory. 

A complete derivation of the equations of relative motion has been presented 
even though only the case of planar motion will be considered in this report. The 
planar and non-planar motions are mathematically and physically distinct problems 
requiring decidedly different techniques for solution. It is felt then that a separate 
consideration of the non-planar motion and it's stability characteristics is justifed. 

In 1960 Clohessy and Wiltshire [ 1 ] published a paper concerned with the 
relative motion of twoclose satellites. Their development concerned itself with 
circular orbits about a spherically symmetric central body. H. S. London [2 ] 
extended the results by including second order terms in various components of 
the relative distance and relative velocity normalized with respect to orbital 
radius and orbital velocity. His work is confined to circular orbits and spherically 
symmetric central bodies. In [3] , J. de Vries made a generalization by allowing 
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the nominal orbit to be slightly eccentric. His work assumes a spnerically symme- 
tric central body. In [ 4 ] , Anthony and Sasaki obt^'n an approximate solution of 
the equations of relative motion incliding linear and quadratic terms in the rela- 
tive distance for the case of nearly circular orbits about a spherical central body. 

In [5], Euler and Shulman made a slight extension of the work of Anthony and 
Sasaki by allowing orbits of arbitrary eccentricity. In [6], Schechter and Cole 
are concerned with the influence of air drag and oblateness. In their study the ori- 
gin is assumed to be located in an infinitely heavy satellite moving in a circular 
orbit at an altitude of 100 n miles. The motion of each of the satellites, obtained 
by an application of a two variable expansion procedure, is referred to this ficticious 
satellite and then, by subtraction, is made relative to each other. 

FORMULATION OF THE EQUATIONS OF RELATIVE MOTION 

The equations of relative motion will be derived from the Lagrangian formu- 
lation where the Lagrangian L is represented by the difference between the kinetic 
energy T and the potential energy U of a spherically symmetric central body. The 
force per unit mass will be denoted by F. 

Choose an inertial coordinate system X Y Z where X is in the plane of the 
equator of the central body and lies along some fixed direction in space, Y is also 
in the plane of the equator ninety degrees ahead of X and Z completes an orthog- 
onal coordinate system whose origin is at the center of mass of the central body 
(see Figure 1). 

Now consider the motion of one of the two satellites as completely known 
relative to this inertial system and define a frame of reference O', whose origin 
is at the center of mass of this satellite, in such a way that x is along the radius 
vector to O', y is in the orbit plane of 0' ninety degrees ahead of x and z is normal 
to the orbit plane of O'. The inclination of the orbit of 0' is denoted by t , it’s 
longitude of the ascending node is 0 and the sum of it’s true anomaly, v, and 
argument of perigee, w, is denoted by u. 

The problem now is to formulate the equations of motion of the satellite M 
with respect to the moving reference frame O'. It is easily seen that 


r 


M 



+ R 


and that 


d r M d V dR 
d t d t + dt 


(1) 


(2) 


2 




where all terms in (2) are referred tc the inertial coordinate system. Taking 
into account that 0' is moving with respect to the fixed inertial system, the deriv 
ative dR/d t can be represented as 



(3) 


where 5, the angular velocity vector of the moving reference frame O', is given 
by 

o>=(i cosu + Qcos u sin l) x 


+ (- £ sin u + ficos u sin t) y 


+ (u + n cos t) z 
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Here the dots indicate derivates with respect time and the symbol * indicates 
a unit vector. Note that 2 is expressed in the moving reference frame. A deri- 
vation of the components of in tne moving coordinate system is presented in 
Appendix 1. 

If the velocity of O' is resolved along the instantaneous directions of the 
moving axes, yielding the components { x f) # , y 0 », z n #}, th<»n the kinetic energy T of 
the M satellite written in the inertial frame, takes the f>rm 


T 




♦ 



Moving 


■f ( a ) X 



(4) 


where m is the mass of the M satellite. Carrying out the necessary operations 
leads to 


t = 5^ x o # + y 0 2 ' * * 0 ^ + ^ (* 2 + y 2 ♦ 

+ ^.{(yaj t - z a > y ) 2 + (z oj k - x u ) g ) 2 ♦ (x o> y - y c^) 2 } 

+ in ly (x - y w + z + y 0 i (y - z ^ + x o) f ) + z q( (z - x w + y oj^)} 
♦ m { x ( - y co g + z <*> y ) + y ( - z co K + x eu g ) + i ( - x w + y ^)> 


(5) 


The first term in this expression is just the kinetic energy of the moving 
origin. The second term is the kinetic energy as would be calculated by an 
observer in the 0' frame who thinks 0' is fixed. The third term is related to 
the moment of inertia about the axis of rotation and may be written simply as 
I w 2 where 


T 


m 

2 


r2 _ (R* 

- " u? 


The fourth term represents the coupling of the motion of 0 # to the motion of M 
and may be rewritten in the form 

m i (R- V) - m (R- A) 
d t 
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where V and A are respectively the velocity and acceleration of 0' . The last term 
may be written in the form 


\ (m (y i - * y)] ♦ cu y [m(zx-xi)] ♦ co g [m (xy-yx)] 

which is recognized as the scalar product oi the angular velocity vector and the 
apparent angular momentum vector, that is, the angular momentum as would be 
calculated by an observer in the 0' frame who thinks O' is fixed. 

The equations of relative motion expressed in generalized coordinates, (q, q), 
take the form 


d_ /BL\ _ BL F 

dt XbaJ " 


( 6 ) 


where L = T - U. Allowing q to take the notation x, y, z respectively yields, from 
the I.agrangian formulation 


x + x Q # + 2 (to z - u> g y) + z ci - y u g + co z Q i - co g y Q 


/ v / x 1 B U P 

+ % ( z % - * & (y % - * %) + - — = F x 

y a y * y m d x x 


(7) 


y + y 0 » + 2 (cu g x <o % z) + x 0 ) f - z co g + eo g x Q « - cu K s ( 


i CO. 


/ v x 1 BU _ 

« < x % - y y ",>+ sr F y 


( 8 ) 


z + z* » + 2 (oj y - co x) + y cL - x co + oj y rt # - co x n 
0 ' x y' ^y y x •'0 y 0 


1 BU 


+ % (y CO - * oj ) - OJ (z O) - x cO + - TT— = F 

y * f x x * m d z 1 


(9) 


Neglecting for a moment the potential U and the components of force it is 
seen that the total acceleration of the M satellite written in the rotating frame 
is made up of three terms: 


(i) acceleration relative to the rotating frame; this leads to terms like 
x , y , * 

(ii) acc< leration due to Coriolis effects; this lead;? to terms like 

2 (<*> y z - co g y), 2 (u g x - oj k z) and 2 (w y - »o y x) 

(iii) acceleration as would be calculated by an observe*: fixed with respect 
to the rotating reference frame; this leads to te^rns like 

V + (co y z q • - cx> t y 0 O + (z o; y - y - x w J + (R- Z) 

V + (^ f * 0 ' - i 0 .) + (x co g - z o; x ) - y co 2 + OJ y (R • Z) 

V + (oj x y Q • - cui y x Q 0 + (y - x ^ y ) - z a; 2 + o; f (R • Zj) 

In order to cast the equations (7) thru (9) into a form which will yield a first 

order theory more readily certain simplifications will be made. To begin <x> x 
will be neglected so as to uncouple the planar motion and the non-planar motion. 

A constraint on a component of the total rotation of the orbital plane of 0' is 
thereby imposed namely, the motion governed by the equation 

£ cos u + fl sin u sin i 

is being neglected. What this means physically is that motion of the orbital plane 
of 0' about the radius vector of O' is being held fixed. Secondly the mutual 
attraction of the two satellites is being neglected. Thirdly it is noted that oj y must 
be indentically zero because the xy-plane must contain, at every instant, the radius 
vector and velocity vector of O' . The preceding comments concerning o> x and a; 
reduce the equations of relative motion to the form 


X + x # . - 2 y W t - y u, - u t y 0 < - x u>\ + ± 

_ 

(10) 

y + y. • + 2 x cu + X u + eu ' x. 1 + y oi? + I 

u z m 

3U- F 
By * 

(11) 

V •• 1 BU r 

z + V + ~ = F z 

u m d z 1 


(12) 
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The formulation of the equations of motion is almost complete. To finish the 
objective certain representations must be used. In particular, the potential U 
takes the form 


u--iL (13) 

m 

where M is the gravitational parameter of the central body and 

r* = (x + x 0 #) 2 + y 2 + z 2 ; 

furthermore the components of force per unit mass due to the oblateness of a 
central body of radius take the form 


F =_u T R 2 

x 2 ^ J 20 e 


F = 


(1-3 sin 2 i sin 2 u) 

i 

0 

(14) 

o2 sin u cos u sin 2 t 
’ • . 4 , 

(15) 

0 


D 2 sin u sin i cos i 

20 

(16) 




Since the motion of 0' is being assumed to obey perturbed tv/o body laws, it is 
clear that the governing equations in the rotating reference irame are 


v - v +4- = {F «v 

* 1 ' 


(17) 


_1 d 

V dt 



= < F yV 


(18) 


where (F x ) 0 * and (F y ) 0 # denote the components of the force on O' due to oblateness. 
In addition, since o> x is being neglected, y 0 ' = x 0 < co g and the angular momentum 
of the perturbed motion will be approximated by that of the two body problem. 
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Making use of (13) thru (18) and the proceeding comments, the system (10) 
thru (12) may be written in the form 


x - 2 y a; + 2 y — - 

1 + 


3 e sin v e co s v 

x u*- 


e cos v 


* 1 + e cos v 


= -r-M J 20 R* -2- (1 - 3 sin 2 t sin 2 u) 


X ^ 9 

*o 


y + 2 x oj m - 2 x 


a>‘ 


e s l n v 


-yaj; 


e cos v 


1 1 i e cos v * 1 + e cos v 


- — M j 5 Q R? — - cos 2 u sin 2 i 
2 2U * 

o 


2 V 




(19) 


( 20 ) 


z + z or 


1 

1 + e cos v 


3 T 
= 2 MJ 


20 


R 


2 sin u sin t cos 
e <■ 


( 21 ) 


These equations represent the complete formulation of the equations of relative 
motion where oblateness of the central body is considered and when the planar 
and non-planar motions are uncoupled, 

SOLUTION OF EQUATIONS GOVERNING PLANAR MOTION 

In this section attention will be devoted to the solution of equations (19) and 
(20). These equations are coupled linear second order differential equations 
containing a small parameter e = 3/2 J 2Q and are amenable to a solution which 
takes the form of a double expansion in eccentricity and in the small parameter. 
As already mentioned the eccentricity is assumed to be small (< 0.01) and 
therefore a first order theory will be presented; consequently consideration will 
be given to the formal summations 

x = x oo CO + e x 01 (t) + € x 1Q (t) (22) 

y (*) = y 00 (O + e y 01 (t) + e y 1Q ( t) (23) 

which are to be substituted into (19) and (20). To derive the equations which are 

satisfied by Xq 0 , x Q1 , . . . , y 10 it is first necessary to expand in terms of e 

and e ; thus 
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( 24 ) 


cj - co + e w + f oj 

* * 0 *i * 


where, in terms of the orbital elements of the rotating origin 0' it is found that 



(25) 


co g = 2 n cos v. (26) 

co g is unknown at this point of the development; it s representation, to be obtained 
later, will serve the purpose of eliminating certain mixed terms which appear 
in the solution of x 10 (t). Furthermore the representation 


( 1 + e cos v)” 1 



■E 


b tl cos k v 
lk 


(27) 


where 




(28) 


will be used in the expansion of the equations of motion. Finally the independent 
variable, time, must be introduce by way of the true anomaly. To accomplish 
this use is made of the Hansen coefficients which are coefficients in the Fourier 
series representation of (r/af exp (i y v) in terms of the mean anomaly, M. 
Lanzano [7] has carried the expansion to second order of eccentricity; the 
expansion to first order is sufficient for this development and is given by 



exp (i yv) = exp i yM - — (a - 2 y) exp (i (y + 


1) M) 


- — (a + 2 y) exp i (y - 1) M 


(29) 


where 


exp i y v = cos y v + i sin y v 


a, y = 0, ± 1, ± 2, • •• • 
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The substitution of (22) thru !29) into (19) and (20) yields the following system 
of equations 


x oo - 2R yoo = 0 


(30) 


X 01 “ 2n ^01 “ 4n y 00 cos n * + 2n 2 y 00 sin n t - n 2 x 00 cos n t =0 (31) 


‘oo 


*10 - 2 "yio -2". Voo = (-7) " 2 x oo( 1 -f sinJ ‘) 


3/ r .\ J -2 

+ 2\T/ " X °° sln 1 


cos (2 W Q + 2n) t 


(32) 


y 0 o + 2 " *oo = 0 


y , 0 + 2S x i 0 + 2 ", x oo =(~) " 2 y 0 <> sin2 ‘ cos (2 w o + 2n) 4 


(33) 


y 01 + 2n Xqi + 4n x 00 cos n t - 2 n 2 x Q() sin n t - n 2 y QQ cos n t = 0 ( 34 ) 


(35) 


where the substitutions M = nt and w = w Q t have been made. 

The circular orbit, spherical central body contribution to the total first order 
solution is governed by equations (30) and (33); the solution is represented as 


x 00 (O = a i sin 


2 nt + a 2 cos 


2nt + 5L 
2n 


(36) 


yoo (O = cos 2 nt - a 2 sin 2 n t + a 3 


(37) 


where a. j = 0, 1, 2, 3 
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are constants which will be evaluated after the x (t) and y (t) solutions are 

10 ' MO 

determined. This solution is almost periodic due to the presence of Cl cos i in 
the expression for the frequency. 

Substitution of (36) and (37) into (34) yields the equation 


n 2 i 8n n 


y oi + 2nx oi = o a 2 sin (2 n + n) t 


3 n 2 - 8 n n 


n 2 + 8 n n 


a 2 sin (2 n - n) t cij cos (2n + n) t (38) 


3 n 2 - 8 n n 


a j cos (2n - n) t + ii 2 a 3 cos nt + n a Q sin nt 


This equation is integrated once and then substituted into (31); carrying out the 
necessary operations leads to the equation for a forced harmonic oscillator 
which x 01 must satisfy, specifically 

x + 4n 2 x Q1 = 4> 01 (t) (39) 

where 4> 01 (t) is given in Appendix II. The solution is given by 

x oi ( t ) = sin n t + (3° 2 l x cos nt + sin (2n + n) t 

+ /3J 2 sin (2n - n) t cos (2n + n) t (^®) 

+ /3°J cos (2n - n) t 

where the j = 1, .... 6 are given in Appendix II. The solution of equation 
(38) now follows easily and is given by 

y Q1 (t) = sin n t + /3® 2 cos n t + sin(2n + n)t 

+ sin(2n-n)t+ /3® 2 cos (2n + n) t (41) 

+ /3° 2 cos (2 n - n) t 
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where the j ■ 1, .... 6 are also to be found in Appendix II. The solutions 
x Q1 (t) and y 0J (t) represent the first order contribution due to the eccentricity 
of the orbit of O'; this contribution is almost periodic. 

Attention will now be turned to equations (32) and (35) which govern the depend- 
ence of the motion on the small parameter. It is noticed that the equations depend 
on the circular orbit spherical central body solutions. It is this dependence 
which leads to mixed terms in the solution; that is, terms of the form 

t sin St, t cos St. 


Specificantly the term 


(t) S2x °<> M sinJ l ) 

leads to terms which oscillate at the same frequency as the solution of the homo- 
geneous equation. Elimination of these mixed terms will be accomplished by 
solving for the small parameter dependence of the angular velocity. 

As before (36) and (37) are substituted into (35) which is then integrated once 
to yield y 10 (t). With an expression for y (t) available it is easy to derive the 
equation which governs the x 10 motion namely, 


x l0 (t) + 4S 2 x,o (t) = <J> 10 (t) (42) 

where <J) i0 (t) is given in Appendix III. It is noticed that(j> 10 (t) contains two terms 
each of which oscillate with a frequency of 2n. Setting each of the coefficients to 
zero to the dependence of oo % on the small parameter; explicitly 


'V 


CO 

x 


1 

8 




(43) 


The solution of (42) is then given by 

x io W = W sin 2 (w Q + n) t + cos 2 (w Q + n) t 

+ sin 2 (w Q + n + n) t + sin 2 (w Q + n - n) t 


( 44 ) 


+ cos 2 (w 0 + n + n) t + /3>« cos 2 (w # + n - n) t + /3* 


10 

x 
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where the/3 J^j = 1, . . . , 7 are to be found in Appendix III. The substitution of 
x 10 (t) into the first integral of (35) yields the solution 


y io = ^!y sin 2 n t + /3>o C os 2n t + /3\° y sin 2 (w () ♦ »0 t 
+ cos 2 (w Q + n) t + /3j° sin 2 (w Q + n ♦ n) t 

(45) 

+ ‘JJ sin 2 (w 0 + n - n) t + ( 3 cos 2 (w Q + n + n) t 
+ /3*° cos 2 (w Q + n - n) t + t 
where the/3j° j * 1, ...» 9 are also found in Appendix III. 


Equations (44) and (45) yield the dependence of the motion upon the small 
parameter, it is noticed that this motion is dependent upon the secular rate of 
change of the true anomaly and that, although the mixed terms have been eliminated, 
a purely secular term remains in the y 10 contribution. The y 10 contribution 
then exhibits an almost periodic behavior superimposed upon a secular drift. 


The complete first order theory for the case of planar motion is represented 
by equations (36), (37), (40), (41), (44) and (45); these representations must be 
combined in the way shown in (22) and (23). The constants of integration 
a. j= 0, 1, 2, 3 may now be easily evaluated; using x (0), x (0), y (0), y (0) to 
denote initial values, it is found that 
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where Aj , Bj , Cj , D , A 2 , B , C 2 , D 2 are functions of the eccentricity and the 
small parameter and are given explicity in Appendix IV. 

CONCLUSION AND COMMENTS 

The first order development of the relative motion of two spacecraft exhibits 
an almost periodic behavior imposed upon a secular drift. The secular contri- 
bution arises when the oblateness of the central body is taken into account; in 
particular, when (32) and (35) are considered. The secular term fails to appear 
in the lower order contribution x Q() (t) and y Qo (t), because the components of the 
derivative of U are restricted to the first order only; this is justified in light of 
the fact that for the application to RAE C & D the ratio of x to x Q « is roughly 
10“ 4 . 


The dependence of the angular velocity on the small parameter disappears 
when the inclination is 54°73. At such an inclination the secular contribution 
also disappears and the motion becomes strictly almost periodic. 

It has been pointed out that the subclass of the total motion which depends 
on the small parameter exhibits a dependence upon the secular rate of change 
of perigee. This dependence in turn leads to a singularity at an inclination 
of 67°8. 

The treatment of the out-of-plane motion will be presented in a subsequent 
report. Furthermore studies are currently being conducted to obtain numerical 
results by comparing the first order analytic theory with the numerical integration 
of equations (7) and (9); these will be reported shortly. This comparison is 
expected to show the effect of the x component of the angular velocity, which, as 
it is recalled, was neglected to uncouple the planar and non-planar motion. 
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APPENDIX I 


DERIVATION OF THE COMPONENTS OF Z IN THE MOVING 

COORDINATE SYSTEM 

Let two position of the moving frame of reference 0' be determined by the 
Euler angles i ,fi, u and i ♦ d • , Q ♦ dQ and u ♦ du. The increment dt corresponds 
to an infinitesimal rotation about the line of nodes ON (see the accompanying 
figure) 


Z 



A 

n 

Figure 2. 


Similarly dD and du correspond to infintesimal rotations about the 0Z and Oz 
respectively. The angular velocity ^ therefore has the components i , 0 , u 
along the axes ON, OZ, and Oz respectively. The decomposition of n , K and k 
into the Oxyz system is easily seen to be 


n = X cos u - y sum 


K = x sin u sin ; 4 y cos u sin 1 ♦ 
k = z 


and therefore 


-* . . • , * 

00 r (t cosn i U sin u sir» l) x 

• • 

+ (- l sin u + flcos u sin 1'iy 
t (u + flcos 1 ) z 
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APPENDIX II 


The forcing function of equation (39) is 

<t> 01 (t) = A.® 1 sin (2 n + n) t + sin (2 n - n) t 

♦ *r cos (2n +n)t + cos (2 n - n) t + cos n t 

+ A.® 1 sin n t 

where 


, oi _ n (4 n + n) (n + 8 n) 

~ “ — a i 

2(2n + n) 1 


id _ n (4n - n) (3n - 8n) 

* 2 ~ = a i 

2 (2 n - n) 1 


x o l _ n (4n + n) (n + 8 n) 

3 ” . . a _ 

2 (2n + n) 2 




n (4 n - n) (3n-8n) 
2 (2 n - n) 




n 


— ( * n - n) a 
2 n 0 


K x - 


n 2 


2 — (n - n) a. 
n 3 


The coefficients appeal ’-g in equation (40) are 


t \ oi _ 6 

r ix 


4 n 2 - n 2 


*Si* 2 


n (4 n - n) 


18 





19 



APPENDIX III 


The forcing function of equation (42) is 


<t> 10 (t) = \j° sin 2 n t + \*° cos 2 n t + \>° sin 2 (w Q 

♦ cos 2 (w Q + n) t + \*° s in 2 (w Q + n + n) t 

+ \*° sin 2 (w Q + n - n) t + \*° cos 2 (w Q + n + 

+ \J° cos 2 (w Q + n - n) t + A.* 0 

where 




"3 (w Q + n) + 5 n 
_ 4 (w Q + n + n) 



3 (w Q + n) + 5 n' 

4 (w Q + n - n) 



"3 (w Q + n) + 5 n “ 
. 4 (w Q + n + n) _ 


+ n) t 


n) t 
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q # 


■©*« 


2 sin 2 t 


3 (w Q + n) - 5 ii 
. 4 (w + n - n) J 


^9° = T (—) n fl - — sin 2 i\ a 

9 4 \a / V 2 0 


The coefficients appearing in equation (44) are 


e\\ = 


q # 




4 [n 2 - (w Q + n) 2 ] 


*i° 


4 [n 2 - (w. + n) 2 ] 


*1! = 


q° 


4 (w + n) [2 n + (w ft + n)] 


/3 10 = 

4x 


C 


4 (w Q + n) [2 n - (w Q + n)] 


p\1- 


q° 


4 (w Q + n) [2 n + (w Q + n)] 


B 10 = 

A 6x 


B 10 = 

h 7x 


K° 


4 (w 0 4- n) [2 n - (w 0 + n)] 

Ss(t) (‘-I' 1 "’ 1 ) 


The coefficients appearing in equation (45) are 


Ct> 




n 
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0 > 


P\° = — a, 

2y n 1 


^° y = 


B 10 = 

4y 


^° y = 


/3 10 = 


/3 10 = 


/3 10 = 


n (w Q + n) M 
4 [n 2 - (w. + n) 2 ] 




2 (w Q + n) 2 - n 2 
4 n [ii 2 - (w Q + n) 2 ] 

n (w n + n)" 1 

4 (w, 


(w 0 + n)-‘ 


R(W o +nr ‘ \ io 1 P.'i*-, • , . “l 

X.1 0 + — l — ) nosin' t 

+ n + n) (2 n + w Q 4 n) 2 W 4 (w + n + n) 2 


n (w Q + n)” 1 


4 (w Q + n - n) (2 n - w Q -n) 


n (w ft + n)" 1 


4 (w Q + n + n) (2 n + w Q + n) 
n (w Q + n)" 1 

4 (w n + n - n) (2 n - w ft - n) 


- j (f) sin J l 


\|® - n J sin J i 


4 (w n + n - n) 


~\2 


\*° ~ ^ (”~) n 2 sin 2 l 


4 (w + n + n) 2 


4 (w. + n - n) 2 


/310 _ 
9y 


"g 

2^ a 0 . 

n 
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APPENDIX IV 


Coefficients appearing in the constants of integration are 


A 1 n 4 n - n 
A. = — - e — + e 

2 n 2n 4 n 2 - n 2 


3 I'W 


n 


16 a ' f» 2 - (w Q + n) 2 

fetr) (R)M { l -l sin2 ) 


sin 2 t 


n , n4n 2 +nn-8n 2 

B = 1 + e 

n 4 n 2 - n 2 


fK\ 2 V 1 

3 (w Q + n) - 5 n 

3 ( w Q + n) + 5 n 

\a / 16 (w Q + n) 

_Cw Q + n -n) [2 n - (w Q +n)] 

(w Q + n + n) (2 n + w Q 


s in 2 l 


l 




n , n 8 i 2 - n n - 16 n 2 j * 

C, = i + e + e \ — 

2 n 


(2 n + n) (2 n - n) 


n 


i. s * p2 1 


L (w Q + n + n) 2 (w Q + n - n) 2 


1 /^e\ n 2 sin 2 l 

' 16 W 


w 0 + n 


3 (w Q + n) n+ 5 n 2 
L (w Q + n + n) 2 (2 n + w Q + n) 


3 (w Q + n) n - 5 n 2 


(w Q + n - n) 2 [2 n - (w + n)]J 


D. = 


n 2 4 n - n n 2 sin 2 i p2 “ 2 ^ w o + 

1 - e — - e| 


A 2 = 2 n + e 


n 4 n 2 - n 2 
2 n 2 


a 


(w f n) 2 4 [n 2 - (w + n) 2 ] 


n 




in 2 i 


w 0 + n 


'3 (w Q + n) + 5 n 3 (w Q + n) - 5 n 
+ 


2 n + w. + n 


2 n - (w + n) J 
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B r 


2 n 2 (n - n) 1 / R e\ 2 n 3 sin 2 t 


+ € — 


4 n 2 - n 2 2 \ a / fj2 - (w n ♦ n) 2 


C 2 = 


_ e n (n - n) e 
4 n 2 - n 2 


3 e\ n 2 s i n 2 i 


8 \ a / -2 


<Ny 

oo. 


n 2 - (w + n) 2 n J 


D 2 = 


r> - 2 n 2 I -'V, lie 

- 2 n - e + e{ - 2 oo m + _ — 


^ \ r* 3 c; i r»2 


n° sin^ i 


n 


1 / R e\ n 3 sin 2 i 

4 8"\T7 


2 \ a / -2 


n 2 - (w Q + n) 2 
3 (w Q + n) + 5 n 


w o +n L( w 0 + n + n) (2 n + w + n) 


3 (w ft + n) - 5 n 


(w Q + n - n) (2 n - (w Q + n) 


► 
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